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Abstract 


A  computationally  fast  algorithm  for  tracking  targets  in  a  sequence  of  scenes 
is  described.  The  algorithm  is  based  on  a  variation  of  template  matching. 
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1.  Introduction 


Among  the  plethora  of  automatic  target  recognition  (ATR)  techniques  de¬ 
veloped  over  the  last  several  decades,  one  of  the  most  basic  is  template 
matching.  This  report  details  a  template-matching  algorithm  developed  to 
track  targets  contained  in  a  sequence  of  forward-looking  infrared  (FLIR) 
images.  The  two  prime  driving  factors  in  algorithm  development  were 
tracking  accuracy  and  computational  burden. 
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2.  Template-Matching  Algorithm 


The  template-matching  concept  is  elementary:  A  rectangular  patch  of  im¬ 
age  pixels  containing  a  representation  of  the  subject  target  is  stepped  across 
a  scene,  and  the  best  match  is  assumed  as  the  most  probable  location  of  the 
target. 

2.1  Scene  Rescaling  and  Preprocessing 

One  can  rescale  the  pixels  of  a  scene  in  a  number  of  ways.  Since  images 
from  the  frame  sequence  for  this  study  are  the  source  of  all  templates  and 
rescaling  involves  some  computational  penalty,  no  frame  rescaling  was 
tried. 

Rather  than  just  using  a  gray-scale  image  as  the  template  basis,  one  can 
preprocess  the  templates  and  frames  to  enhance  various  image  attributes 
and  deemphasize  others.  Although  any  preprocessing  scheme  can  increase 
computational  requirements,  its  use  can  be  justified  where  significant  per¬ 
formance  gains  can  be  achieved.  As  a  first  try,  a  gradient  model,  based 
on  the  Sobel  operator,  was  implemented.  Initial  results  indicated  that  the 
gradient-based  template  performed  so  much  worse  than  the  gray-scale- 
based  template  that  all  further  investigations  into  preprocessing  were 
ended. 

2.2  Background  Correlation 

The  primary  source  of  frame-to-frame  target  location  variability  for  the  test 
scene  set  is  random  frame  dither.  To  estimate  target  location  in  any  given 
frame  with  respect  to  background  in  the  previous  frame,  one  can  use  back¬ 
ground  correlation.  For  the  test  scene  set,  the  maximum  frame-to-frame 
variation  of  target  location  (based  on  ground-truth  data)  in  either  the  hori¬ 
zontal  or  vertical  directions  is  10  pixels. 

To  compensate  for  frame  dither,  I  used  a  rectangular  template  centered  in 
any  given  frame,  with  a  15-pixel  border.  For  the  128-  x  128-pixel  images, 
the  corresponding  template  encompasses  approximately  10,000  pixels — an 
exceedingly  large  template.  To  vary  such  a  template  over  a  conservatively 
chosen  range  of  15  pixels  is,  from  a  computational  point  of  view,  imprac¬ 
tical.  I  found  that  such  a  template  could  be  down  sampled  by  selection 
of  every  fourth  pixel  vertically  and  horizontally  (resulting  in  a  background 
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template  of  approximately  600  pixels)  without  adversely  affecting  perform¬ 
ance  accuracy.  The  template  can  also  be  shifted  in  steps  of  two  pixels  over 
the  15-pixel  range — resulting  in  a  further  factor  of  four  reduction  in  back¬ 
ground  correlation  time.  These  reductions  in  computational  requirements 
allow  template  matching  to  become  a  practical  approach  to  aligning  se¬ 
quential  frames. 

2.3  Scene  Search  Region 

With  the  given  or  computed  location  of  a  target  in  any  frame,  by  using 
background  correlation,  one  can  estimate  the  target's  location  for  the  next 
frame.  This  would  seem  to  permit  a  target  search  near  the  region  of  its 
estimated  location.  Such  a  localized  search  will  reduce  algorithm  execution 
time  and  may  also  improve  performance. 

To  search  over  increasingly  smaller  regions  of  a  frame  while  minimizing 
the  probability  of  the  target  falling  outside  the  search  region,  I  included  a 
frame-to-frame  variable  search  region  capability  into  the  algorithm.  To  im¬ 
plement  this  capability,  I  chose  a  nominal  search  region  size  for  the  data  set. 
With  the  actual  or  computed  target  location  in  any  frame  (from  the  back¬ 
ground  correlator),  the  target's  location  in  the  next  frame  was  estimated.  If 
after  the  most  probable  location  of  the  target  on  the  second  frame  (within 
the  given  search  region)  was  computed  from  the  target  template  and  if  the 
estimated  and  calculated  locations  agreed,  the  search  region  size  for  the 
next  frame  continued  unchanged.  If  this  distance  between  the  estimated 
and  calculated  locations  grew,  then  the  search  region  for  the  subsequent 
frame  could  be  enlarged  (to  a  maximum  search  area  four  times  the  nomi¬ 
nal  search  area).  If  the  calculated  location  of  a  target  approached  too  closely 
to  the  border  of  the  search  region,  the  probability  of  error  for  that  frame  was 
assumed  to  have  significantly  increased.  Because  so  many  frames  were  in 
the  sequence  of  images,  not  all  were  required  to  track  a  target;  these  high 
probability  error  frames  were  discarded,  and  the  search  continued  to  the 
next  frame.  I  found  that  the  cumulative  effect  of  these  procedures  improved 
performance  while  reducing  execution  time. 

2.4  Template  Updating  and  Weighting 

If  you  were  to  track  any  target  through  a  sequence  of  frames,  its  appearance 
would  progressively  change.  This  change  would  have  the  obvious  effect  of 
degrading  the  performance  of  the  template-matching  algorithm.  To  com¬ 
pensate  for  this  degradation,  at  least  in  part,  the  pixel  values  of  the  target 
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template  are  progressively  updated  as  the  algorithm  passes  down  the 
sequence  of  frames.  Equation  (1)  describes  how  these  values  are  updated: 

T&+ i(n)  =  (Pi T&{n)  +  FA(n))/(Pi  +  1)  ,  (1) 


where 


Pi  =  user-defined  constant, 

T(n )  =  the  value  of  pixel  n  of  template,  and 
F(n)  =  the  value  of  pixel  n  for  the  calculated  location  of  target  A 
in  frame  K. 


It  is  reasonable  to  assume  that  not  all  pixels  of  a  template  contribute  equally 
to  a  template-matching  algorithm's  performance.  For  one  to  determine 
whether  this  unequal  contribution  is  random  from  frame  to  frame  or  is  (in 
some  manner)  systematic,  template  weighting  was  introduced.  If  this  con¬ 
tribution  is  not  highly  random,  then  template  weighting  should  improve 
performance.  To  implement  the  weighting,  simply  multiply  the  results  of 
the  template-frame  match  for  each  template  pixel  by  the  weight  term.  This 
weight  is  computed  as  (one  should  note  the  similarity  between  equations 
(1)  and  (2)) 


Wtf+i («)  =  (foWfon)  +  aA(n)aA(n))/(p2  +  1)  ,  (2) 

where 

aA(n)  =  (S  -  sA(n))/S  , 

S  =  sA(n)  ,  and 

n 

sA(n)  =  (TA(n)-FA(n))(TA(n)-FA(n)). 

I  selected  this  weight  term  after  experimentation  with  several  different 
forms  of  these  equations. 


3.  Algorithm  Performance 


The  performance  of  the  template-matching  algorithm  was  optimized  and 
then  tested  with  the  test  scene  set.  The  L1816  test  set  has  ground-truth  val¬ 
ues  for  all  targets  for  frames  53  through  288.  Parameter  optimization  and 
testing  were  done  with  this  frame  sequence.  This  sequence  can  be  divided 
into  two  regions:  (1)  the  beginning,  where  comparatively  few  pixels  are  on 
target  and  the  frame-to-frame  variation  is  minimal  and  (2)  the  end,  where 
both  the  target  pixel  count  and  frame  variations  are  much  larger.  I  designed 
an  optimum  model  for  each  region  (model  1  and  model  2).  An  optimum 
model  is  one  that  both  maximizes  performance  and  minimizes  algorithm 
execution  time.  I  also  tried  to  keep  as  many  parameters  of  each  region's 
model  the  same  as  possible  (see  table  1). 

The  parameters  are  defined  in  terms  of  the  variable  names  used  in  the  code 
listed  in  the  appendix.  Where  parameters  were  previously  given  alternate 
designations,  these  designations  are  given  in  parentheses  (see  table  1). 


3.1  Model  1 

This  model  is  the  parameter  set  designed  to  perform  optimally  for  the 
frames  at  the  beginning  of  the  frame  sequence.  The  performance  of  model  1 
is  demonstrated  by  defining  a  frame  sequence,  selecting  the  first  frame  as 


Table  1 .  Algorithm 
parameters.  Variable 
names  are  those  of  code 

Model  1 

Model  2 

sample 

= 

1 

sample 

= 

2 

sample 

= 

2 

listed  in  appendix  (see 

bgnd 

= 

4 

bgnd 

= 

4 

bgnd 

= 

4 

function  main). 

stepa 

= 

1 

stepa 

= 

1 

stepa 

= 

1 

step 

= 

2 

step 

= 

2 

step 

= 

2 

reset 

= 

1 

reset 

_ 

1 

reset 

= 

1 

edge 

= 

13 

edge 

= 

20 

edge 

= 

22 

error 

= 

1 

error 

= 

1 

error 

= 

1 

errora 

= 

1 

errora 

- 

1 

errora 

- 

1 

mask 

= 

9 

mask 

= 

11 

mask 

= 

12.5 

weight 

= 

1 

weight 

= 

1 

weight 

= 

1 

wt  (j32) 

= 

6 

wt  (02) 

= 

6 

Wt  (02  ) 

= 

6 

time  (00 

= 

4 

time  (00 

= 

4 

time  (0i ) 

= 

4 

suppress 

= 

0 

suppress 

= 

0 

suppress 

= 

0 
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the  source  of  all  target  templates,  and  testing  its  performance  on  the  subse¬ 
quent  frames.  Two  targets  are  on  each  frame:  the  M60  and  tnk.  To  increase 
the  statistical  significance  of  the  results,  I  repeated  model  l's  demonstration 
of  performance  for  a  series  of  40  frame  sets.  Each  frame  set  was  created  by 
incrementing  the  previous  frame  sequence  (and,  hence,  also  the  target  tem¬ 
plate  source)  by  one  frame.  For  instance,  if  the  first  frame  set  encompassed 
frames  55  through  115  (with  the  template  source  as  55),  then  the  second  set 
would  encompass  56  through  116  (with  the  template  source  as  56),  and  so 
on  for  40  sets.  The  results  for  the  40  frame  sets  were  then  either  averaged 
or  summed. 

To  optimize  the  template-matching  algorithm's  parameters,  I  used  a  20- 
frame  set  with  each  spanning  100  frames.  The  initial  frame  set  encom¬ 
passed  frames  55  through  155.  To  test  the  performance  of  the  optimized 
model,  I  used  the  aforementioned  sequence  of  40-frame  sets.  These  sets 
spanned  varying  frame  sequence  lengths.  Because  both  optimization  and 
test  data  sets  were  drawn  from  the  same  database,  I  sought  to  introduce 
differences  into  the  two  sets.  Besides  the  differences  in  frame  set  counts, 
varying  frame  sequence  starting  points  (hence  different  target  templates) 
were  used.  Nevertheless,  one  should  recognize  the  overwhelming  similari¬ 
ties  between  the  optimization  and  test  data  sets  and  be  very  circumspect  in 
interpreting  the  results.  As  more  data  sets  are  brought  on-line,  this  problem 
should  diminish. 

Results  for  model  1  are  presented  for  two  different  target  templates.  The 
first  is  for  a  template  down  sampled  by  selection  of  every  fourth  pixel  (table 
2)  and  the  second  is  for  no  down  sampling  (table  3).  The  limited  number 
of  pixels  on  target  for  the  data  sets  used  with  model  1  complicates  down 
sampling  and  maintaining  performance  accuracy. 


Table  2.  Performance  of 
down-sampled  model  1. 
Code  parameters  listed  in 
table  1. 


A1 

B 2 

C3 

D 

4 

E5 

M60 

tnk 

M60 

tnk 

both 

55 

120 

33 

16.7 

26.5 

27 

13 

12 

55 

100 

25 

12.9 

16.1 

28 

16 

15 

55 

80 

17 

8.9 

8.9 

29 

23 

21 

55 

60 

9 

4.2 

4.6 

29 

27 

24 

75 

40 

5 

7.8 

2.8 

25 

33 

25 

95 

20 

1 

2.3 

0.7 

30 

37 

30 

'Beginning  frame  number  of  first  frame  sequence. 
2Frame  sequence  length. 

3Average  number  of  skipped  frames  per  sequence. 
4  Average  number  of  bad  target  hits  per  sequence. 
5Total  number  of  error-free  passes  through  se¬ 
quence  for  40-frame  sequences. 
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Table  3.  Performance  of 
model  1  with  no  target 
down  sampling.  Code 
parameters  listed  in 
table  1. 


3.2  Model  2 


A1 

B 2 

C3 

D 

4 

E5 

M60 

tnk 

M60 

tnk 

both 

55 

120 

3 

2.6 

29.8 

38 

2 

2 

55 

100 

2 

2.6 

12.9 

38 

11 

11 

55 

80 

1 

1.9 

3.3 

38 

28 

28 

55 

60 

1 

0.9 

0.9 

38 

38 

38 

75 

40 

3 

0.6 

4.4 

37 

32 

32 

95 

20 

1 

0.1 

1.5 

39 

34 

34 

beginning  frame  number  of  first  frame  sequence. 
2Frame  sequence  length. 

3  Average  number  of  skipped  frames  per  sequence. 

4  Average  number  of  bad  target  hits  per  sequence. 
5Total  number  of  error-free  passes  through  se¬ 
quence  for  40-frame  sequences. 


Note  that  for  models  1  and  2,  the  target  template  is  significantly  larger  than 
the  target.  The  inclusion  of  substantial  background  in  the  template  while 
adversely  affecting  computer  execution  time  significantly  improved  accu¬ 
racy.  Despite  the  large  templates,  where  a  match  was  achieved,  the  differ¬ 
ence  between  estimated  and  ground-truth  location  of  targets  on  average 
was  typically  one  or  two  pixels  and  with  the  worst-case  fit  being  rarely 
more  than  six  pixels. 


Model  2  was  optimized  for  the  frames  near  the  end  of  the  data  set  sequence. 
Its  performance  is  given  in  table  4.  I  used  20  frame  sets  to  optimize  the 
parameters  of  this  model.  The  frame  sequence  for  the  first  of  these  frame 
sets  was  185  through  245.  Template  updating  has  a  significant  effect  on 
model  performance.  This  is  shown  in  table  5. 

The  test  configuration  was  that  of  row  4  of  table  4  with  parameter  time  (/?i) 
the  variable.  Note  that  template  weighting  minimally  affects  performance 
and  can  probably  be  deleted  from  the  model. 
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Table  4.  Performance  of 
down-sampled  model  2. 
Code  parameters  listed 
in  table  1. 


Table  5.  Effect  of 
template-updating 
parameters  (/3i)  on 
performance  of  model  2. 


A1 

B2 

C3 

D 

4 

E5 

M60 

tnk 

M60 

tnk 

both 

65 

180 

— 

20.0 

10.6 

7 

8 

7 

85 

160 

79 

10.9 

6.6 

27 

19 

18 

105 

140 

28 

1.8 

3.7 

38 

27 

27 

125 

120 

22 

0.0 

0.8 

39 

37 

37 

145 

100 

22 

0.0 

0.0 

40 

40 

40 

165 

80 

22 

0.0 

0.0 

40 

40 

40 

185 

60 

22 

0.0 

0.0 

40 

40 

40 

beginning  frame  number  of  first  frame  sequence. 
2Frame  sequence  length. 

3  Average  number  of  skipped  frames  per  sequence. 

4  Average  number  of  bad  target  hits  per  sequence. 
5Total  number  of  error-free  passes  through  se¬ 
quence  for  40-frame  sequences. 


Time1 

c2 

D3 

E4 

M60 

tnk 

M60 

tnk 

both 

10,000 

78 

16.4 

8.8 

12 

17 

12 

20 

68 

19.7 

13.4 

11 

14 

10 

10 

29 

2.2 

15.9 

32 

9 

9 

8 

21 

3.5 

12.2 

27 

16 

16 

6 

20 

2.1 

3.8 

31 

30 

30 

4 

22 

0.0 

0.8 

39 

37 

37 

3 

22 

0.3 

0.7 

37 

36 

36 

4See  appendix  (or  j3\  of  eq  (1))  for  definition  of 
time. 

2Average  number  of  skipped  frames  per  sequence. 
3  Average  number  of  bad  target  hits  per  sequence. 
4Total  number  of  error-free  passes  through 
sequence  for  40-frame  sequences. 
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4.  Discussion 


It  is  feasible  to  develop  both  a  reasonably  accurate  and  a  computation¬ 
ally  efficient  algorithm  for  target  tracking  based  on  template  matching.  The 
frame-to-frame  redundancy  in  a  sequence  of  images  is  an  attribute  that  is 
exploited  by  the  present  version  of  template  matching  to  maximize  track¬ 
ing  accuracy.  As  has  been  demonstrated,  much  that  is  potentially  computa¬ 
tionally  inefficient  about  a  template-matching  scheme  can  be  circumvented 
through  careful  algorithm  design.  As  is  too  often  the  case  in  the  target 
recognition  field,  a  definitive  statement  cannot  be  made  about  absolute  al¬ 
gorithm  accuracy.  Defining  an  absolute  performance  standard  requires  not 
only  a  more  extensive  set  of  image  sequences  than  was  available  for  this 
study  but  also  several  alternative  tracking  algorithms  for  comparison. 
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Appendix.  Target-Tracking  Algorithm  Source  Code 


The  following  is  a  C  source  code  listing  of  the  template-matching  algorithm 
of  this  report. 

/*  template. c  template  matching  code  for  the  SDF  data  sets  */ 

iinclude  <stdio.h> 

Sinclude  <stdlib.h> 
iinclude  <math.h> 
include  <string.h> 

#def ine  size  300 
#def ine  fsize  128 
Idefine  t count  2 

int  seq,  start,  stop,  gtx [tcount] [size] ,  gty [tcount] [size] ,  alt,  counter,  dot; 

int  sample,  row,  col,  rate,  edge,  tplate [tcount] [50]  [50] ,  error,  errora,  ed; 

int  x [tcount] [size],  y [tcount] [size],  tgtx[size],  tgty[size],  sing,  dcount; 

int  xshift,  yshift,  weight,  time,  bgnd,  step,  reset,  bad [tcount],  stepa,  skip; 

int  wrong [tcount] ,  recover[tcount] ,  suppress,  wt ,  box,  skip_seq [200] ; 

int  inc,  pass,  sta,  sto,  kount [t count], efree [tcount] [500]; 

float  mse[tcount] [fsize] [fsize],  back  [tcount] ,  target [tcount] ,  mask,  dist; 

float  frame! [2] [fsize] [fsize],  mask! [tcount] [size] ,  w[tcount] [50] [50] ,  bigd; 

float  bad.d [tcount]; 

char  isource[7],  image  [12]; 

double  frame [fsize] [fsize]; 

unsigned  char  imtfsize]  [fsize] ; 


read.ground.truthO 

{ 

int  na,  nb,  nc,  ka,  xx,  yy,  big; 

char  stringa[12],  stringb [5] ,  stringc [12] ,  temp [81]; 

FILE  nn.filel,  *in_file2; 

for(nb=0;  nb<tcount;  nb++){ 
for(na=0;  na<size;  na++){ 
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gtx [nb] [na] =-500 ; 
gtytnb]  [na]=-500; 
}} 


if (seq==l ){ 

in_f  ilel=fopen("/net/ragu/export/data/sdf/motion/L18i6/L1816m60 .gt2" ,nr") 
inJile2=fopen(7net/ragu/export/data/sdf/motion/L1816/L1816tnk.gt2Vr") 
/*  presently  have  gt  values  in  above  files  for  frames  053  thru  288  */ 
for(na=0;  na<1000;  na++){ 
for(nc=0;  nc<2;  nc++){ 
if(nc==0){ 

fscanf (in_f ilel , "'/.s  */,s  Xd  XdIT,&stringa?  fotringb,  fox,  &yy); 
if(feof(injilel))goto  jumpa; 

/*  following  line  stores  enough  characters  to  later  find  image  source  */ 
if(na==0)for(nb=0;  nb<7;  nb++) isource [nb] =stringa [nb] ; 
for(nb=0;  nb<5;  nb++)stringb[nb]=stringa[nb+7] ; 
ka=atoi (stringb) ; 

/*  xx  and  yy  contain  x  and  y  ground  truth  locations...  note  corrections  to 
the  values  stored  in  the  ground  truth  tables  7 
gtx[nc]  [ka]  =xx; 
gty [nc] [ka]=fsize-yy; 
fgets(temp,  80,  in.filel); 

} 

if (nc==l){ 

fscanf (in^file2, "Xs  '/,s  Xd  Xd'r,&stringc,  fotringb,  fox,  ftyy); 
if (strcmp(stringa,stringc) !=0){ 
printf ("PROBLEM  WITH  GROUND  TRUTH  DATA\n"); 
exit(l); 

} 

gtx[nc]  [ka]  =xx ; 

gty [nc] [ka] =fsize-yy; 

fgets(temp,  80,  in.file2); 

} 

} 

} 

jumpa: 

fclose(injfilel); 
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fclose(in_file2); 


/*  find  maximum  frame-to-f rame  axis  shift  for  gt  values:  */ 
big=0; 

for(na=0;  na<t count;  na++){ 
for(nb=53;  nb<288;  nb++){ 
ka=nb+i; 

if  Cabs(gtx[na]  [nb]  -gtx  [na]  [ka])>big)big=abs(gtx[na]  [nb]  -gtx  [na]  [ka] ) ; 
if  (abs(gty  [na]  [nb]  -gty  [na]  [ka]  )>big)big=abs(gty  [na]  [nb]  -gty  [na]  [ka]); 
}} 

printf ("largest  f rame-to-frame  gt  axis  shift=  ’/,d\n\big) ; 

} 

} 


get _f rame (na) 
int  na; 

{ 

int  nb,  nc,  ka=0,  kb=0; 
double  big,  small; 

char  temp[81],  string[100],  stringa[4]  ; 
FILE  *in_file,  *out_file; 


if (seq==l){ 

if  (na<10)sprintf  (strings,  "000,/,d",na); 

if (na>9  kk  na<100)sprintf (strings, "00’/,d",na); 

if (na>99  kk  na<1000)sprintf (strings, "0'/,d\na); 

if(na>999)sprintf(stringa,",/,d">na); 

strcpy(image,isource); 

strcat (image, strings); 

/*  image  []  contains  frame  na  file  name  */ 

} 

sprintf  (string/'/net/ragu/export/data/sdf /motion/Locaas/7,  s.rl.bin",  image); 
inJile=fopen(string,"rb") ; 
if (na==start){ 
big=-l.E10; 
small=l.E10; 
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} 

for(nb=0;  nb<3;  nb+^)fgets(temp ,  80,  in.file); 

for(nb=0;  nb<fsize;  nb++){ 

for(nc=0;  nc<fsize;  nc++){ 

fread(&im[nb] [nc] , sizeof (unsigned  char) , 1 , in.file) ; 

frame  [nb]  [nc]= (double  )im[nb]  [nc] ; 

framel [0] [nb] [nc] =framel [1] [nb]  [nc] ; 

framel  [1]  [nb]  [nc]=frame[nb]  [nc] ; 

if (na=-start){ 

if (frame [nb] [nc] <small) small=f rame [nb] [nc] ; 
if (frame  [nb]  [nc] >big)big=f rame [nb] [nc] ; 

} 

}} 

fclose(in.file); 

h  can  use  following  block  to  generate  a  single  MATLAB  image  */ 

/* 

if (na==start){ 

out.file=fopen(,,scene.dat",,lw"); 
for(nb=0;  nb<f size ;  nb++){ 
for(nc=0;  nc<f size ;  nc++){ 

fprintf (out.f ile , "Id  " , (int ) (100 . * (big-frame [nb] [nc] )/(big-small) ) ) ; 
} 

fprintf  (out  .file, 'An1'); 

} 

fclose(out.file); 

} 

*/ 

} 


down_sample_frame(ka) 
int  ka; 

{ 

int  na,  nb,  ma,  mb; 
ma=-l; 

for(na=0;  na<f size;  na++){ 
if(na'/,sample==0){ 
ma+=l; 


mb=-l; 

} 

for(nb=0;  nb<f size;  nb++){ 
if  (nay,sample==0  &&  nb'/,sample==0){ 
mb+=l; 

frame [ma] [mb] =f rame [na] [nb] ; 

} 

}} 

/*  row  and  col  contain  the  row  and  column  count  for  the  down  sampled  frame:  */ 


alt  ernat  e .do wn.s  ampl e .frame ( ka) 
int  ka; 

{ 

int  na,  nb,  ma,  mb; 


ma=-l; 

for(na=0;  na<fsize;  na++){ 
if  (nay,sample==0){ 
ma+=l; 


} 

for(nb=0;  nb<fsize;  nb++){ 
if(nay,sample==0  &&  nb'/. s amp  1  e== 0 )  { 
mb+=l; 

if (na>0  &&  na<f size-1  kk  nb>0  k&  nb<f size-1) 
frame  [ma]  [mb]  =0 . 4*f rame  [na]  [nb]  +0 .  l*f rame  [na-1]  [nb]  +0 .  l*f rame  [na+1]  [nb] + 
0. 1+frame  [na]  [nb-l]+0.1*frame[na]  [nb+l]+0.05^frame  [na-1]  [nb-l]  + 

Q.05*f  rame  [na-1]  [nb+l]+0.05*f  rame  [na+1]  [nb-l]+0.05*frame  [na+1]  [nb+1] ; 
else 

frame  [ma]  [mb]  =frame  [na]  [nb] ; 

} 

}} 

/*  row  and  col  contain  the  row  and  column  count  for  the  down  sampled  frame:  */ 
if(ka==start+l){ 
row=ma+l;- 
col-mb+1 ; 


} 
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templatesJroirLgt(ka) 
int  ka; 

{ 

int  na,  nb,  nc,  ma,  mb,  ia,  ib; 
for(nc=0;  nc<tcount;  nc++){ 

/*  following  block  allows  the  first  frame  thru  template  matching  to  be  masked  */ 
if (ka==start  il  reset==l){ 
x[nc]  [ka]=gtx[nc]  [ka] /sample; 
y  [nc]  [ka]=gty[nc]  [ka] /sample; 

} 

ma=gty[nc] [ka] /sample-edge/2; 
if (ma<0)ma=0; 

if((ma+edge)>col)ma=col-edge; 
mb=gtx[nc] [ka] /sample-edge/2; 
if(mb<0)mb=0; 

if((mb+edge)>row)mb=row-edge; 

ia=-l; 

for(na=ma;  na<ma+edge;  na++){ 

ib=-l; 

ia«l; 

for(nb=mb;  nb<mb+edge;  nb++){ 
ib+=l; 

/*  tplatefnc] [ia] [ib] =( int ) (0 . 5*(f ramel [0] [na] [nb]+framel[l] [na] [nb]));  */ 
tplate [nc] [ia] [ib]=frame [na] [nb] ; 

}} 

} 

} 


do_t emplate.mat ching (nc , naa) 
int  nc,  naa; 

{ 

int  na,  nb,  nd,  ka,  kb,  ma,  mb,  ja,  jb,  jc,  jd,  xx,  yy; 
float  aa,  m,  temp,  low,  ax,  ay,  a; 
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low=l.E+40; 

m=2.0*mask*inaskl[nc]  [naa] ; 
for(na=0;  na<row;  na++){ 
for(nb=0;  nb<col;  nb++){ 
mse  [nc]  [na]  [nb]  =-1.0; 

}} 

j  a=( int ) (y  [nc]  [naa-1] +yshif t/sample-row/m) ; 
jb=(int)(x[nc] [naa-1] +xshift/s ample -col/m) ; 
jc=(int)(y[nc]  [naa-1] +yshif t/sample+row/m) ; 
jd=(int)(x[nc] [naa-1] +xshif t/ sample+col/m) ; 


for(na=ja;  na<=jc;  na+=stepa){ 

if (na<edge/2+l  II  na>=row-l -edge /2) goto  jumpa; 

for(nb=jb;  nb<=jd;  nb+=stepa){ 

if(nb<edge/2+l  1 1  nb>=col-l-edge/2)goto  jumpb; 

temp=0 .0; 

aa=0.0; 

roa=-l; 

f or (ka=na-edge/2 ;  ka<na-edge/2+edge;  ka++){ 

ma+=l ; 

mb=-l; 

for(kb=nb-edge/2;  kb<nb-edge/2+edge;  kb++){ 
mb+=l; 

if (error==0){ 

temp+=w  [nc]  [ma]  [mb]  *fabs(tplate  [nc]  [ma]  [mb]  -frame  [ka]  [kb] ) ; 
aa+=l; 

} 

if (error==l){ 

temp+=¥[nc]  [ma]  [mb]  *  (tplate  [nc]  [ma]  [mb] -frame  [ka]  [kb])* 

(tplate  [nc]  [ma]  [mb]  -frame  [ka]  [kb] ) ; 

aa+=l; 

} 


} 

} 

mse  [nc]  [na]  [nb]=temp/aa; 
if  (mse  [nc]  [na]  [nb]<low){ 
low=mse  [nc]  [na]  [nb] ; 


/*  make  next  search  box  a  function  of  location  of  target  in  present  box:  */ 
ax=fabs((x[nc]  [naa-1] +xshift/sample-xx)/(col/(inask*maskl  [nc]  [naa]))); 
ay=fabs((y[nc] [naa~l] +yshif t/sample-yy )/(row/(mask+maskl [nc] [naa] ) ) ) ; 
if (ax>ay) 
a=ax; 

else 

a=ay; 

h  printf(  "counter:  '/,d  a['/,d]=  7,.3f  ",  counter  ,nc  ,a  [nc] ) ;  */ 

if(a>0.45  M  suppress==0){ 
skip=l; 

/*  printf ( "SKIPPING  FRAME  #  7,d\n",naa);  */ 
goto  jump; 

} 

/*  printf ("\n");  */ 
maskl  [nc]  [naa+l]=(1.0-a); 
x[nc]  [naa]=xx; 
y  [nc]  [naa]=yy; 
jump: 


predict.target.locationO 

{ 

int  na,  nb,  ma,  mb,  ka,  kb,  kc; 
float  temp,  tempa; 

h  use  background  correlation  to  find  frame  to  frame  dither 
ka  is  maximum  anticipated  axis  shift 
kb  is  template  matching  step  size 
kc  relates  to  pixel  count  in  background  template  */ 
temp=l.E+40; 


ka=15; 

kb=step; 

kc=bgnd; 

for(na=-ka;  na<=ka;  na+=kb){ 
for(nb=-ka;  nb<=ka;  nb+=kb){ 
tempa=0 .0; 

for(ma=ka;  ma<fsize-ka;  ma+=kc){ 
f or(mb=ka;  mb<fsize-ka;  mb+=kc){ 
if(errora==0) 

tempa+=fabs (framel  [1]  [ma]  [mb]  -framel  [0]  [ma+na]  [mb+nb] ) ; 
if (errora==l) 

tempa+=(framel  [1]  [ma]  [mb] -framel  [0]  [ma+na]  [mb+nb])* 
(framel  [1]  [ma]  [mb]  -framel  [0]  [ma+na]  [mb+nb] ) ; 

}} 

if(tempa<temp){ 

temp=tempa; 

xshift=nb; 

yshift=na; 


}} 


} 


detected_target_weighting(naa) 
int  naa; 

{ 

int  na,  nb,  nc,  ma,  mb; 

float  temp[tcount]  [50]  [50] ,  tempa[tcount] ,  a; 

for(nc=0;  nc<t count;  nc++){ 

tempa  [nc]=0.0; 

ma=-l; 

for(na=y [nc]  [naa] -edge/2;  na<y [nc]  [naa] -edge/2+edge ;  na++){ 

ma+=l; 

mb=-l; 

for(nb=x[nc] [naa] -edge/2;  nb<x[nc]  [naa] -edge/2+edge ;  nb++){ 
mb+=l ; 

temp[nc]  [ma]  [mb]=l  .0; 

if (na>=0  M  na<row  M  nb>=0  &&  nb<col) 
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temp  [nc]  fma]  [mb] =(tplate  [nc]  [ma]  [mb]  -frame  [na]  [nb])* 
(tplate [nc] [ma] [mb] -frame [na] [nb] ) ; 
tempa [nc] +=temp [nc] [ma] [mb] ; 

}} 

} 

for(nc=0;  nc<tcount ;  nc++){ 
for(na=0;  na<edge;  na++){ 
for(nb=0;  nb<edge;  nb++){ 
a=(tempa[nc]-temp[nc] [na] [nb] ) /tempa [nc] ; 
w[nc]  [na]  [nb]=(wt*w[nc]  [na]  [nb]+a*a)/(wt+l.G); 

}}} 

/*  adjust  below  term  to  appropriate  value  to  print  weight  set  */ 
if  (naa7,4000==0)  { 

for(nc=0;  nc<tcount;  nc++){ 
for(na=0;  na<edge;  na++){ 
for(nb=0;  nb<edge;  nb++){ 
printf 3f  " ,w[nc]  [na]  [nb] ) ; 

} 

printf ("\n"); 

} 

printf("\n\n"); 

} 

} 

} 


generate^averaged.template(naa) 
int  naa; 

{ 

int  na,  nb,  nc,  ma,  mb,  t; 


t=time; 

/*  if ( counter?, rate==l)t=l;  */ 

for(nc=0;  nc<tcount;  nc++){ 
ma=-l; 

for(na=y[nc] [naa] -edge/2 ;  na<y[nc] [naa] -edge/2+edge ;  na++){ 
ma+=l ; 


mb=-l; 

for(nb=x[nc] [naa] -edge/2;  nb<x[nc] [naa] -edge/2+edge ;  nb++){ 


if (na>=0  &&  na<ro¥  kk  nb>=0  kk  nbccol) 

tplate[nc]  [ma]  [mb] = ( t  *tplat  e  [n  c]  [ma]  [mb]+frame  [na]  [nb])/(t+1.0); 

}} 

} 

} 


generate_performance_statistics(ka,kb) 
int  ka,  kb; 

{ 

int  nc; 

float  temp,  er; 


dot=0 ; 

for(nc=0;  nc<t count;  nc++){ 
if  (counter— 1 )  wrong [nc] =0 ; 
er=sqrt( 

(1.0*gtx[nc]  [ka]/sample-x[nc]  [ka])*(1.0*gtx[nc]  [ka]/sample-x[nc]  [ka])+ 
(1.0*gty[nc]  [ka]/sample-y[nc]  [ka] )*(1  .O^gty [nc]  [ka] /sample-y [nc]  [ka])) 
if (er>0.5*edge){ 
efree[nc]  [kb] =1 ; 
bad[nc]+=l; 
bad_d[nc]+=er; 
dot=l; 

if (wrong [nc] =0) kount [nc] +=1 ; ' 
wrong [nc]=l; 

/*  printf ("bad  hit  on  */,d  at  counter  '/.d  error/edge=  X.3f  edge=  ’/,d\n" 
,nc, counter, er/edge, edge);  * / 


else{ 

if ( (counter- l)Xrate ! =0  kk  wrong [nc]==l){ 
re cover [nc]+=l; 
wrong [nc]=0; 

} 

wrong [nc] =0 ; 

temp=sqrt((l .  0*gtx[nc]  [ka]/ sample -x[nc]  [ka] )+ 
(1.0*gtx[nc]  [ka] /sample-x  [nc]  [ka])+ 


( 1 . 0*gty [nc]  [ka] / sample-y [nc] [ka] ) * ( 1 . 0*gty [nc] [ka] / sanple-y [nc]  [ka] ) ) 
dist+=temp; 

if (temp>bigd)bigd=temp ; 
dcount+=l; 

} 

} 

} 

make  ATLAB_raovi e (naa) 
int  naa; 

{ 

int  na,  nb,  nc,  ma,  ka,  kb,  temp[fsize] [fsize] ; 
float  big,  small,  m; 
char  string [100] ; 

FILE  *out_file; 

big=-l.E10; 

small=l.E10; 

for(na=0;  na<row;  na++){ 
for(nb=0;  nb<col;  nb++){ 
if(frajne[na]  [nb]  >big)big=f  rame  [na]  [nb] ; 
if (frame[na] [nb] <small)small=frame [na] [nb] ; 

}} 

for(na=0;  na<row;  na++){ 
for(nb=0;  nb<col;  nb++){ 

temp [na]  [nb] =100. *( (f rame [na] [nb]-small)/(big-small)) ; 

}} 

/*  box  gt  location  of  targets:  */ 
for(nc=0;  nc<tcount;  nc++){ 
m=2 . *mask*maskl [nc] [naa]; 
ma=100; 

if(nc==l)ma=20; 

f or(na=gtx[nc] [naa] /sample -edge /2~1 ; 

na<gtx[nc] [naa] / sample-edge/2+edge+l ;  na++){ 
temp[gty[nc]  [naa]  /saraple+edge/2+1]  [na]=ma; 
temp[gty [nc] [naa]/sample+edge/2-edge-l] [na]=ma; 


} 

f or(na=gty [nc] [naa] / sample-edge/2-1 ; 

na<gty[nc]  [naa] /sample-edge/2+edge+l ;  na++){ 
temp  [na]  [gtxfnc]  [naa]  /  sample+edge/2+1]  =ma; 
temp [na] [gtx[nc] [naa] / sample+edge/2-edge-l] =ma; 

} 

/*  cross  hairs  on  best  template  matching  fit:  */ 
for(nb=0;  nb<6;  nb++){ 

if  ((y[nc]  [naa]  -edge/2-nb)>=0)temp[y  [nc]  [naa] -edge/2-nb]  [x[nc]  [naa]]=ma; 
if  C(y  [nc]  [naa]+edge/2+nb)<row)temp[y[nc]  [naa] +edge/2+nb]  [x[nc]  [naa]]=ma; 
if  ((x[nc]  [naa]-edge/2-nb)>=0)temp  [y  [nc]  [naa]]  [x[nc]  [naa] -edge/2-nb]  =ma; 
if((x[nc]  [naa]+edge/2+nb)<col)temp[y[nc]  [naa]  ]  [x  [nc]  [naa]+edge/2+nb]=ma; 

} 

/*  dotted  line  around  search  box  for  each  target  based  on  value  of  mask  */ 
if(mask>=0.01  kk  box==0){ 
for(na=x[nc] [naa-l]+xshift/sample-col/m; 

na<x [nc] [naa-1] +xshift/sample+col/ni;  na+=2){ 
if (na>=0  kk  na<col){ 

if ( (y [nc]  [naa- 1] +yshif t /sample+row/ m) <row) 

temp[(int)(y[nc] [naa- 1] +y shift / s  ample+row/m) ] [na]=ma; 
if ( (y [nc] [naa- 1] +yshif t /sample-row/m) >=0) 

temp  [(int) (y [nc] [naa-1] +yshift/sample-row/m)] [na]=ma; 

} 

} 

f or (na=y [nc] [naa-1] +yshif t /sample-row/m; 

na<y[nc] [naa- 1] +y shift/ sample+row/m ;  na+=2){ 
if (na>=0  kk  na<row){ 

if ( (x[nc] [naa- 1] +xshif t /sample* col/m) <col) 

temp[na] [(int)(x[nc] [naa-1] +xshif t/sample+col/m)]  =ma; 
if ( (x [nc] [naa- 1] *xshif t /sample-col/m) >=0) 

temp[na]  [(int)(x[nc] [naa-1] +xshift/sample-col/m)]=ma; 

} 

} 

} 

/*  if  bad  hit  on  frame,  place  large  dot  in  corner  */ 
if (dot==l){ 

for(ka=0;  ka<8;  ka++){ 


for(kb=0;  kb<8;  kb++){ 
temp [ka] [kb]=ma; 

}} 

} 

} 

sprintf  (string ,  "scene  .dat7,d" ,  counter) ; 
out_file=fopen(string,V); 
for(na=0;  na<row;  na++){ 
for(nb=0;  nb<col;  nb++){ 
fprintf (out_f ile/’Xd  n  ,temp[na]  [nb] ) ; 

} 

fprintf (out_f ile , u\n") ; 

} 

fclose(out.file); 

} 


get.singleJATLAB_iniage() 

{ 

/*  creates  side-by-side  down  sampled  frame  and  template  matching  response 
with  targets  (from  gt)  boxed  */ 

int  na,  nb,  nc,  nd,  temp [tcount+1] [f size] [f size] ; 
float  big,  small; 

FILE  *out_file; 

for(nc=0;  nc<tcount;  nc++){ 

big=-l.E10; 

small=l.E10; 

for(na=0;  na<row;  na++){ 
for(nb=0;  nb<col;  nb++){ 
if (mse [nc] [na] [nb] >big)big=mse [nc] [na] [nb] ; 
if (mse [nc] [na]  [nb] <small  M  mse [nc] [na]  [nb] >=0 . 0) small=mse  [nc] [na] [nb] ; 
}} 

for(na=0;  na<row;  na++){ 
for(nb=0;  nb<col;  nb++){ 
if (mse[nc] [na] [nb]>=0.0) 

temp [nc] [na]  [nb]=100 . *( (mse [nc] [na]  [nb] -small) /(big-small)); 


else 

temp[nc]  [na]  [nb]=50.; 

}} 

} 


bip-l.ElO; 
small5 1. El 0; 

for(na=0;  na<row;  na++){ 
for(nb=0;  nb<col;  nb++){ 
if (frame [na] [nb] >big)big=f rame [na] [nb] ; 
if (frame [na] [nb] < small )small=f rame [na]  [nb] ; 

}} 

for(na=0;  na<row;  na++){ 
for(nb=0;  nb<col;  nb++){ 

temp[tcount] [na] [nb] =100 . * ( (frame [na] [nb] -small)/ (big-small)) 

}} 

for(nd=0;  nd<tcount+l ;  nd++){ 
for(nc=0;  nc<t count;  nc++){ 
f or(na=gtx [nc] [sing] /sample-edge/2-1 ; 

na<gtx[nc] [sing]/sample-edge/2+edge+l;  na++){ 
temp [nd] [gty [nc] [sing]/ sample +edge/2+l] [na] =100 ; 
temp[nd] [gty [nc] [sing]/sample+edge/2-edge-l] [na] =100 ; 

} 

for(na=gty[nc] [sing] /sample-edge/2-1; 

na<gty [nc] [sing] /sample-edge/2+edge+l ;  na++){ 
temp[nd] [na] [gtx [nc] [sing] /sample+edge/2+l] =100; 
temp[nd] [na] [gtx[nc] [sing] /sample+edge/2-edge-l] =100; 


} 

} 

out.f ile=fopen("scene .dat" ,Mff") ; 

for(na=0;  na<row;  na++){ 

for(nc=0;  nc<tcount+l;  nc++){ 

for(nb=0;  nb<col;  nb++){ 

fprintf (out Jile/Xd  ",temp[nc] [na] [nb]); 

} 

} 


} 


fprintf(out_file,"\n"); 

} 


mainO 

{ 

int  na,  nb,  nc,  nd,  ka,  kb[tcount],  kc,  ja,  jb,  ma,  quit,  total; 
FILE  *out_f ile ; 

printf ("image  data  set  sequence?  l=1816\n"); 
scanf  ("*/,d,T  ,&seq) ; 

printf ("start  and  end  sequence  of  frames  at  frames  number :\n") ; 

scanf (M'/,d  JSd"  ,ftstart  ,ftstop) ; 

sta=start; 

sto=stop; 

printf ("increment  size  between  passes  thru  data  set?\n") ; 
scanf  Old", Sine); 
pass=l; 
if (inc>0){ 

printf ("number  of  passes?\n"); 
scanf  ("7,d"  ,&pass) ; 

} 

printf ("down  sample  frames  by  selecting:\n"); 

printf ("every  pixel  (1)  every  second  (2)  every  third  (3)...?\n"); 

scanf  ("7td"  ,&sample) ; 

if (sample>l){ 

printf ("select  down  sample  mode:  A=  1  B=  2\n"); 
scanf  Old", ftalt); 

} 

printf ("corresponding  sample  rate  for  background  template :\n"); 
scanf  ("7,d"  ,&bgnd) ; 

printf ("foreground  pixel  stepping  rate :\n") ; 
scanf  Old", festepa); 

printf ("background  pixel  stepping  rate : \n") ; 
scanf  Old", ftstep) ; 


printf ("reset  target  to  gt  value  for  each  new  template?  no=0  yes=l\n") ; 
scanf  ("Kd\ftreset); 

printf  ("square  templates:  edge  length  in  pixels  before  down  sampling?\n"); 

scanf  (’’Xd"  ,&edge) ; 

ed=edge; 

edge/=sample; 

printf ("in  down  sampled  pixels:  \nu); 

printf  ("form  of  template  match?  error=0  squared  error=l  \n11); 

scanf  Old",  terror); 

printf ("and  for  background:  error=0  squared  error=l\n"); 
scanf(’"/,d",&errora); 

printf  ("mask  frame  away  from  predicted  target  location\n"); 
printf ("no=  0.5  yes=  divide  axis  lengths  by : \n" ) ; 
scanf("Kf",fcnask); 

printf ("include  template  pixel  weighting?  no=0  yes=l\n"); 
scanf  (,f7,d"  ,&weight) ; 
if (weight==l){ 

printf ("weighting  time  constant  (int)?\nM); 
scanf  ("7,d"  ,&wt) ; 

} 

printf ("time  constant  for  template  averaging  (int):\nu); 
scanf  ("Xd\fttime); 

printf ("1st  frame  always  selected  in  the  sequence  as  a  template  source\n"); 
printf ("template  selection  rate?  (e.g.  5=  every  fifth  frame)\n"); 
scanf  ("7.d"  ,&rate) ; 

printf ("suppress  frame  skipping?  no=0  yes=l\n"); 
scanf ("Xd" ,&suppress) ; 

printf  ("plot  MATLAB  single  frame  with  corresponding  template  response\n"); 
printf  ("and  gt  locations?  no=  0  yes=  frame  number  or,  alternatively, \n"); 
printf ("make  (down  sampled)  frames  for  MATLAB  movie:  input  -l\n"); 
scanf  (l,Xd"  ,&sing) ; 
if(sing==-l){ 

printf ("remove  search  box  from  movie  frames?  no=0  yes=l\n"); 
scanf ("Xd" ,&box) ; 

} 

read_ground_truth() ; 
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}} 

for(nc=0;  nc<t count;  nc++){ 
for(na=0;  na<50;  na++){ 
for(nb=0;  nb<50;  nb++){ 
w[nc]  [na]  [nb] =1.0; 

}}} 

if (sample==l){ 
row=fsize; 
col=fsize; 

} 

if (sample>l){ 

row=f size/sample+l; 
col=row; 

} 

get.frame (start); 

if  ( sample>l )  down.samplejrame  ( start ) ; 

templates_from„gt (start) ; 

count er=0 ; 

na=start; 

skip=0; 

jb=0; 

for(ja=start+l;  ja<=stop;  ja++){ 
j  umpa : 

if(sldp==l)jb+=l; 

na+=l; 

if(na>=stop){ 
skip.seq[jb]+=l; 
goto  juinpb; 

} 

get.frame (na); 

predict  .target.locationQ ; 

if (sample>l  alt==l)down.sample_frame(na); 

if (sample>l  &&  alt==2)alternate.down_sampleJraine(na); 

skip=0 ; 

for(nd=0;  nd<tcount;  nd++){ 
do.template .matching (nd,na) ; 
if (skip==l){ 

for(nb=0;  nb<t count;  nb++){ 


maskl [nb] [na+l]=0.5; 
x  [nb]  [na]=x[nb]  [na-1] ; 
y[nb]  [na]  =y  [nb]  [na-1]; 

} 

for(nb=0;  nb<f size ;  nb++){ 
for(nc=0;  nc<f size ;  nc++){ 
framel  [1]  [nb]  [nc]=framel  [0]  [nb]  [nc] ; 

}} 

goto  jumpa; 

} 

} 

total+=l; 

counter+=l; 

skip_seq[jb]+=l; 

jb=0; 

if(veight==l)detectedJ;argetjreighting(na); 
generate_averaged_template(na) ; 
generate_perforinance_statistics(na,ma); 
if (sing>0  &&  na==sing)get.single_HATLAB.iinage(); 
if (sing==-l)make_MATLAB_movie(na) ; 
if  (counter*/,  rate==0)templates_from.gt(na); 
if(na>=stop)goto  jumpb; 

} 

jumpb : 


for(na=0;  na<tcount;  na++)kb[na]=0; 
for(na=0;  na<t count;  na++){ 
for(ma=0;  ma<pass;  ma++){ 
if(efree[na]  [ma]==0)kb[na]+=l; 

}} 

ka=0; 

for(ma=0;  ma<pass;  ma++){ 
kc=0; 

for(na=0;  na<tcount;  na++){ 
if (efree [na] [ma] !=0)kc=l; 

} 

if (kc==0)ka+=l; 
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} 


if (seq==l  M  stop==288) 

printf ( "WARNING . . .may  have  exceeded  limit  with  parameter:  stop\n"); 
pnntf ("following  results  are  average  values  per  pass  thru  sequence^11); 
printf  ("number  of  error  free  passes  thru  sequence=  7.d\n",ka); 
printf  ("total  frame  count  per  pass=  7, .  if\n" ,  (f  loat)total/pass) ; 
for(nc=0;  nc<tcount;  nc++){ 

printf  ("target  ‘/,d  number  of  error  free  passes  thru  seq.=  MV,  nc,  kb  Inc]); 
printf  ("target  M  bad  hit  count=  l/(.lf\n"Jnc,(float)bad[nc]/pass); 
printf  ("target  M  avg  bad  hit  pixel  distance  from  gt=  7,.lf\nu,nc, 
bad.d[nc]/(bad[nc]+0.001)); 

printf  ("target  '/,d  recovered  bad  sequence  count=  7..2f\n",nc, 
(float)recover[nc]/pass); 

printf  (’’target  */,d  bad  sequence  count=  7,.2f\n",nc, 

(float)konnt [nc]/pass) ; 

} 

printf  ("avg.  detected  target  gt  distance  in  down  sampled  frame=  '/,.2f\n", 
dist/dcount); 

printf  ("large  st  detected  target  -  gt  distance=  '/,.2f\n",bigd); 
printf  ("total  target  hit  count  per  pass=  7,d\n",dcount/pass); 
printf  ("consecutive  frame  skip  sequence  length  vs  total  population\n"); 
for(nc=0;  nc<200;  nc++) if ( skip.seq [nc] >0 ) printf ( "7.d=7,d  ",nc,skip_seq[nc]); 
printf ("\n\n"); 

out_file=fopen( "template . dat" , "w" ) ; 

f printf (out .file, "test  performance  of  template. c  (17Sept98  version)\n"); 
fprintf (out .file,"  pass  7.d\n  inc  7ld\n">pass,inc); 
f  printf  (out  .file,"  seq  7,d\n  start  (initial)  */,d\n  stop  7*d\n  sample  M\n  bgnd  */,d 
step  7,d\n  reset  */,d\n  edge(original)  M\n  error  M\n  weight  '/,d\n  time  M 
mask  7, .  3f  \n" ,  seq ,  st  a,  s  to ,  sample ,  bgnd ,  st  ep ,  re  set ,  ed ,  error ,  weight ,  t  ime ,  mask) ; 
if(sample>l)fprintf(out.file,"  alt  7#d\n" ,alt) ; 
fprintf  (out  .file,"  errora  M\n  rate  7.d\n",errora,rate); 
fprintf  (out  .file,"  stepa  7*d\n  suppress  7id\n",stepa,  suppress); 
if(weight==l) 

fprintf  (out  .file,"  wt  7.d\n",wt); 

fprintf (out.f ile , M\n  total  frame  count  per  pass=  7»d\n", total/pass); 
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fprintf (out.file, "following  data  is  average  value  per  pass  thru  sequence\n"); 
fprintf (out .file, "number  of  error  free  passes  thru  sequence1 7,d\n’1  ,ka) ; 
for(nc=0;  nc<tcount;  nc++){ 

fprintf (out.file, 11  target  '/,d  number  of  error  free  passes  thru  seq.=  Xd\n" , 
nc  ,kb[nc] ) ; 

fprintf  (out.file,"  target  '/,d  bad  hit  count=  ,Mf\n",nc, 

(float)bad [nc] /pass) ; 

fprintf  (out.file, 11  target  ‘/,d  recovered  bad  sequence  count=  X  2f\nlf  ,nc, 
(float)recover[nc]/pass); 

fprintf  (out.file, 11  target  */d  bad  sequence  count=  l2f\n",nc, 

(float)kount [nc]/pass); 

} 

fprintf (out.file, 

"avg.  detected  target  gt  distance  in  down  sampled  frame=  ,/,.2f\n",dist/dcount); 
fprintf  (out.file ,  "largest  detected  target  -  gt  distance=  '/,.2f\n\bigd); 
fclose(out.file); 
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